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Coherent manipulation of a Majorana qubit by a mechanical resonator 
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We propose a hybrid system composed of a Majorana qubit and a nanomechanical resonator, 
implemented by a spin-orbit-coupled superconducting nanowire, using a set of static and oscillating 
ferromagnetic gates. The ferromagnetic gates induce Majorana bound states in the nanowire, which 
hybridize and constitute a Majorana qubit. Due to the oscillation of one of these gates, the Majorana 
qubit can be coherently rotated. By tuning the gate voltage to modulate the local spin-orbit 
coupling, it is possible to reach the resonance of the qubit-oscillator system for relatively strong 
couplings. 
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I. INTRODUCTION 

Majorana bound states (MBSjP^ are recently attract¬ 
ing increasing interest both theoretically and experimen¬ 
tally. These have been predicted to exist in artificial 
structures, such as nanowires with spin- orbi t coupling 
(SOC) in proximity to a superconductor,^^ ferromag¬ 
netic atom chains on top of a superconductor,!^ topo¬ 
logical insulator/superconductor hybrid structures,^®^! 
and superconducting circuits.^i^ Recently, possible s ig- 
natures of MBS have been reported in nanowires 
atom chains,!^ and topological insulator/superconductor 
structures.li^ The MBS attract considerable atten¬ 
tion partly due to their hypothetical non-abelian 
anyonic statistics, which might allow the realiza¬ 
tion of top ologi cally-protected quantum information 
manipulation.l2iiSH2ni In parallel to the on going search of 
some unambiguous confirmation of MBS,!2il^ there are 
also numerous theoretical studies on how to efhciently 
exploit these MBS. 

One promisiM application of MBS is to construct Ma¬ 
jorana qubits.ES It has been suggested that Majorana 
qubits might be robust against local perturbati ons and 
are hence promising to store quantum information.^I^SM] 
(Note that Majorana qubits are not totally protected 
from decoherence, as studied in, e.g., Refs. 
Furthermore, Majorana qubits could be rota ted by 
topologically-protected braiding operations.^^^ISll There¬ 
fore, among various realizations of qubits,Majo¬ 
rana qubits are considered to be promising candidates for 
building blocks of quantum information processors. The 
braiding operations alone are insufficient to realize a uni¬ 
versal quantum gate based on a Majorana qubit.^^ For 
the implementation of arbitrary qubit rotations, other 
non-topological operations are required. Several schemes 
of su ch non -topological op eratio ns assisted by, e.g ., phase 
gates,!23111 quantum dots,l^2E2l flux qubits, SRS or mi¬ 
crowave cavities,!^ have been proposed in the literature. 

Nanomechanical resonatord^H could also be used to 
study non-topological operations of a Majorana qubit. 
For example, quite recently, Kovalev et have pro¬ 


posed to rotate a Majorana qubit by a magnetic can¬ 
tilever. Indeed, nanomechanical resonators have been 
utilized to couple to a wide range of quantum systems, 
includiM electric circuits,!^ optomechanical devices,!^ 
atoms,Cooper-pair boxes,!^ spin qubits,!^ or mi¬ 
crowave cavities.lSil With the assistance of nanomechani¬ 
cal resonators, it is possible to perform important appli¬ 
cations such as quantum manipulations, quantum mea¬ 
surements, as well as efficient sensing. These applications 
exploit the advantages of nanomechanical resonators, 
e.g., their large quality factors (10^-10®), high natural 
frequencies (MHz-GHz), as well as the feasibility of reach¬ 
ing the quantum ground states by cooling methods.!^^^!^ 
Recently, nanomechanical resonators have also been ex¬ 
ploited to measure or manipulate the MBS.I^^lt^ Never¬ 
theless, the study of hybrid system^ coupling nanome¬ 
chanical resonators to Majorana qubits is quite limited. 
This work aims to contribute to this field. In this paper, 
we propose another Majorana qubit-nanomechanical res¬ 
onatorhybrid system in the framework of the spin-boson 
model,!^ based on a semiconductor nanowire in proxim¬ 
ity to an s-wave superconductor. We show that a strong 
coupling between a nanomechanical resonator and a Ma¬ 
jorana qubit can be achieved, allowing an efficient trans¬ 
fer of quantum information between these two quantum 
systems. Further, with braiding operations, it should be 
possible to realize a universal quantum gate based on a 
Majorana qubit. 

This paper is organized as follows. First, we describe 
the Majorana qubit and its coupling to a nanomechani¬ 
cal resonator. Afterwards, we numerically study the cou¬ 
pling strength and the resonance condition of the hybrid 
system. Then, we solve the qubit-phonon dynamics and 
achieve a coherent control of the Majorana qubit. Fi¬ 
nally, we summarize our results. 

II. MODEL AND HAMILTONIAN 

As illustrated in Fig. we consider a semiconductor 
nanowire with a Rashba SOC of strength on the sur¬ 
face of an s-wave superconductor with a superconducting 
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gap A. Three ferromagnetic gates, FMl, FM 2 and FM 3 , 
are placed on top of and along the nanowire. Among 
these gates, FMl and FM 3 are static while FM 2 is free 
to harmonically oscillate along the nanowire (with a mass 
M and an oscillation frequency ujq). The gates FMl 
and FM 3 are sufficiently long (of the order of 1-10 /im) 
while FM 2 in between is relatively short (of the order of 
100 nm). These ferromagnetic gates induce a local Zee- 
man splitting in the nanowire. For simplicity, we take 
the Zeeman splitting under the three gates to be identi¬ 
cal, with a magnitude of Bq. An electric voltage V can 
be applied on the gates to modulate the Rashba SOC lo¬ 
cally, e.g., from to Oy. In our study, we consider the 
case with Bq > + where /j, is the chemical poten¬ 

tial in the nanowire. Therefore in the nanowire, the parts 
subject to the Zeeman splitting (under the three ferro¬ 
magnetic gates) are in the topological (T) region. The 
remaining parts, without t he Ze eman splitting, are in the 
non-topological (N) region .EESl As a result, the nanowire 
has an N-T-N-T-N-T-N domain structure where the 
three T domains are under the gates. At the six bound¬ 
aries between the N and T domains in the nanowire, 
MBS arise. As the two outer MBS are far apart, only the 
four inner ones, schematically labeled as 71-74 in Fig. 
are coupled due to hybridization arising from their small 
separatior|S 2 H 63 | and are hence relevant to our considera¬ 
tion. 

To lowest order, the hybrid system constructed above 
can be described by the Hamiltonian 

H = Hm + H,sc, (1) 

where the mutual coupling Hamiltonian of the MBS^SHSI 


Hm = *ffn[^12(i)]7l72 + *54(^23)7273 + *5n[^34(<)]7374, 

( 2 ) 

and the nanomechanical oscillator Hamiltonian 

Hose = + 2'^^o^o(0' (3) 

The coupling strengths gn,t depend on the domain 
lengths kj. Due to the oscillation of the gate FM 2 , 
= li2 + Xo{t) and l^it) = ^34 — Xo{t) are time 
dependent. Here, Xo{t) stands for the displacement of 
the gate FM 2 from its balance position, which is much 
smaller than the static domain lengths 34- Therefore, 
to first order in xq, one has 

Hm =*[5n(^?2) + 2;o(i)5n(^i2)]7i72 + i9t{l23)l2ls 

+ *[5n(^34) - a:^o(05^(^34)]7374- (4) 

The four MBS, satisfying {71,7^} = S ^, c an be used 
to construct a Majorana qubit as follows.^ES At first we 
define two Dirac fermion operator^ = (71 + 
and = (72 -I- i^3)/V^. The Hilbert space of Hm can 
then be spanned by states |n^, n^), with the fermion oc¬ 
cupation numbers = c|c^ and = c|c^. Due to the 
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FIG. 1 : (Color online) (a) Schematic diagram of the proposed 
Majorana qubit-nanomechanical resonator hybrid system. A 
semiconductor nanowire is placed on the surface of an s-wave 
superconductor. Three ferromagnetic gates are on top of the 
nanowire, of which FMl and FM 3 are sufficiently long (of the 
order of 1-10 pm) and static, while FM 2 is relatively short (of 
the order of 100 nm) and free to oscillate as a harmonic oscil¬ 
lator. The ferromagnetic gates induce a local Zeeman split¬ 
ting Bq in the underlying nanowire, and can also be used to 
modulate the local Rashba SOC strength by applying an elec¬ 
tric voltage V. (b) Wave amplitude Ivltp of the four coupled 
MBS at a static state in an InSb nanowire with the set-up 
shown in (a). The red dashed and red solid curves respec¬ 
tively correspond to the wave amplitudes of the lowest two 
eigenstates (close to the zero energy). These two states con¬ 
stitute the Majorana qubit. The dotted curve with the scale 
on the right-hand side of the frame indicates the profile of the 
inhomogeneous Zeeman splitting along the nanowire. In the 
calculation I12 = I34 = 150 nm, I23 = 400 nm, Bq = 1 meV, 
and A = 0.5 meV. The gate voltage V is zero and the Rashba 
SOC is homogeneous along the nanowire, with a strength 
= 20 meV nm. 



conservation of fermion parity, the states {|0,1), |1, 0)} 
and || 0,0), |1 ,1)| form two decoupled (odd and even) 

sectors.EEffil yyg 

assume that there is no high-energy 
excitation (e.g., no Cooper-pair breaking in the super¬ 
conducting substrate) and restrict our study to the odd 
sector with -I- = 1 . For convenience, we define 

pseudo-spins | f) = |1)0) and | |) = |0,1), and use them 
as the two logical states of the Majorana qubit 
In this pseudo-spin space, 17172 = —*7374 = —cry and 
*7273 = —cTz- The nanomechanical oscillator is quan- 
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tized in the Fock space {|n)} with \n) 


^| 0 ), where 


® “ V (^0 + annihilation operator of 

phonons. Consequently, in the space {| ti) ® l^)}j the 
hybrid system can be simply described by the spin-boson 
Hamiltonian,!^ 


HeS =--(Jz - Say + g{a^ + a)ay + nhujQ, (5) 

with the constant omitted. Here e = 2gt{l23), S = 
9 nili 2 ) - Qnilli), and g = -ioK(^i 2 ) + ffn(^ 34 )]. where 
xo = [h/{ 2 MLOo)Y/‘^ is the zero-point motion of the oscil¬ 
lator. 


III. HYBRIDIZATION OF MAJORANA BOUND 
STATES 

In this section, we study the MBS and their mu¬ 
tual coupling. In the static state, the inhomoge¬ 
neous nanowire can be described by a tight-binding 
model. Using the Bogoliubov-de Gennes basis ^j — 
(/iT^ fih where fjy stands for the fermion op¬ 

erator of a spin -?7 {g =t, i) electron on the j-th lattice 
site, the particle-hole Hamiltonian read^ 

Hbdg = ^ + H.c.)], (6) 
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where 

hj — (2to /4 )soT 2 -t- ^sqTx -j- BjSyTQj (7) 

ij = toSoTz + iajSzTz. ( 8 ) 

In the above Hamiltonian, the Pauli matrices Tx,y,z act 
on the particle-hole space and Sx,y,z act on the real 
spin space. The spin-diagonal hopping energy is to = 
h?/{2m*a^), and the spin-off-diagonal hopping energy is 
aj = of/{2a). Here Bj and are the on-site Zeeman 
splitting and Rashba SOC, respectively, m* is the effec¬ 
tive electron mass, and a is the lattice spacing in the 
discretized tight-binding model. In the T {N) domains 
Bj = Bq {Bj = 0) and = ay (a?- = a{/). When the 
gate voltage V is zero, ay = ao • 

Here, to lowest order, we follow Refs. 05] and (ni to 
investigate the coupling strength {gt) approximately 
in an isolated T-N-T {N-T-N) three-domain structure. 
In such a simplified model, the inner N {T) domain has 
a finite length, while the outer two T {N) domains are 
assumed to be infinitely long. By numerically diagonal¬ 
izing this three-domain system, one can obtain the en¬ 
ergy splitting of the two MBS localized at the two T/N 
boundaries. This energy splitting is precisely caused by 
the coupling of the MBS. With and gt known numeri¬ 
cally, the Majorana qubit can be well described by e and 
(5, and the qubit-phonon coupling g can be obtained also 
from g/ [refer to Eq. Moreover, by exactly diagonal¬ 
izing the Hamiltonian of the genuine N-T-N-T-N-T-N 


domain structure as shown in Fig. [^a), one can obtain 
the hybrid four MBS under consideration. 

In this work, we consider an InSb quantum wir^i^with 
an effective electron mass m* = 0.015 nie, a Rashba SOC 
a^ = 20 meV nm, and a large Landau factor g^ ~ 50. 
We choose the superconducting gap A = 0.5 meV, the 
local Zeeman splitting Bq = 1 meV, the chemical poten¬ 
tial 13 = 0, and the lattice constant a = 10 nm. The total 
lattice site number is chosen as 1000 for the numerical 
convergence. In Fig. we show the dependence of the 
Majorana coupling strength {gt) on the length of the 
N (T) domain In {k), as well as the derivative g/ ver¬ 
sus In- Further, as an example, in Fig. we present 
the wave amplitude j'l'p of the four hybrid MBS, when 
^12 = ^34 = 150 nm and I 23 = 400 nm. In Fig. [^b), the 
red dashed and red solid curves stand for the wave am¬ 
plitudes of the lowest two eigenstates (close to the zero 
energy) in the static inhomogeneous nanowire. The state 
corresponding to the red solid (dashed) curve is mainly 
contributed by the 72 and 73 (71 and 74) MBS. Here, 
these two states form the Majorana qubit. 



C(i) (nm) 


FIG. 2: (Color online) Majorana coupling strength g„ {gt) 
versus l„ {It), the length of the inner N {T) domain between 
the two outer T {N) domains. The derivative g/ versus In is 
also shown, with the scale on the right hand side of the frame. 
The necessary parameters for the calculation are specified in 
the main text. 


IV. QUBIT-PHONON COUPLING AND 
RESONANCE 

We now look into the qubit-phonon coupling and the 
resonance condition. We assume that the nanomechan¬ 
ical oscillator FM2 has a mass M = 10“^® Kg and an 
oscillation frequency Wq = 5 MHz. With these parame¬ 
ters, the zero-point motion of the oscillator is calculated 
to be Xo = 0.1 pm. We consider the symmetric case 
with li 2 = Z 34 = 150 nm, and hence we have (5 = 0 and 
g = 0.2 MHz in Eq. ([^. The longitudinal length I 23 of 
the FM2 gate is chosen as 400 nm, such that the Rabi 
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resonance condition e « —cuq can be easily satisfied, e.g., 
by further subtly adjusting the gate voltage V which con¬ 
trols the local Rashba SOC strength ay. In Fig. we 
present the variation of e as well as g versus ay. It is 
shown that when slightly adjusting V, and hence a^, 
the resonance point s « —ojq can be reached while the 
qubit-phonon coupling g remains almost invariant. This 
qubit-phonon coupling is relatively strong, in view of the 
long lifetime of the Majorana qubit and the high qual¬ 
ity factor of the nanomechanical oscillator. In principle, 
the qubit-phonon coupling can be stronger when the do¬ 
main length li 2 (as well as ^ 34 ) becomes smaller (refer 
to Fig. [^. However, if the two edge modes 71 and 72 
(as well as 73 and 74 ) are too close and hence their hy¬ 
bridization becomes quite strong, the model Hamiltonian 
([^ describing four distinguishable MBS might fail. 





FIG. 3: (Color online) Qnbit splitting e and the qubit-phonon 
coupling g (the scale is on the right-hand side of the frame) 
versus ay, the Rashba SOC strength in the topological (T) 
domains modulated by the gate voltage V. 


V. QUBIT-PHONON DYNAMICS 

Here we study the dynamics of the qubit-phonon hy¬ 
brid system. To achieve this, we m ake use of the Python- 
based Qutip software packag^SHMI solve the Lindblad 
master equation, 

m = - ^[i?eff,p(0] + 1 E 

k 

+ [Lkp{t),Ll]Y (9) 

In this equation, p is the density matrix of the qubit- 
phonon system, and are the Lindblad operators ac¬ 
counting for the dissipation of the hybrid system due to 
its coupling to the environment. The relaxation of the 
Majorana qubit is taken into account by Li = ^JXjTxa 
while the dissipation of the nanomechanical resonator is 


included by = \/(n + Ijwo/Qa and L 3 = a/ uloqIQo^. 
Here n = [exp(Iia;o/A:BT) — 1]“^ is the thermal phonon 
number in equilibrium with the environmental tempera¬ 
ture r, Q is the quality factor of the nanomechanical os¬ 
cillator, and Ti is the usual relaxation time of the qubit. 
By solving the master equation, one can obtain the time 
evolution of the qubit and phonon occupations. 

In our model, the temperature T is set as 10 mK and 
hence the thermal phonon number n is as large as 258. 
Therefore, an additional cooling of the oscillatoi!^^^^ is 
required, e.g., as also applied in a proposed nanomechan¬ 
ical resonator-nitrogen-vacancy center hy brid system.l^ 
We assume that after side-band cooling^SHU the phonons 
thermally occupy the lowest several quantum states with 
a small phonon number, e.g., n = 0.3. The initial state of 
the Majorana qubit is set as | f)) implying that a single 
electron is splitted into the 71 and 74 Majorana fermions. 
Experimentally, this initial state might be realized when 
only the FMl and FM3 gates are in proximity to the 
nanowire before inserting the middle FM2 gate. The re¬ 
laxation time of the Majorana qubit depends on the con¬ 
crete set-up and environment. Following Refs. and E71 
we typically set Ti around 100 ps. 

In Fig.|^ we plot the time evolution of the occupations 
of the qubit and phonons respectively, with different val¬ 
ues of Ti and Q. As indicated by the figure, quantum 
information can be effectively transferred back and forth 
between the Majorana qubit and the nanomechanical res¬ 
onator. During this process, the single electron in the 
nanowire alternatively occupies (back and forth) the pair 
of MBS: 7 i and 74 , or 72 and 73 . Inversely, this quan¬ 
tum information transfer can also modulate the motion of 
the oscillator, e.g., the oscillation amplitude. In fact, as 
the nanomechanical resonator is near its quantum ground 
state, the oscillation amplitude (xg), which might be ob¬ 
servable, is almost linearly related to the phonon number. 
This is because (xg) oc ((a'^ -|- a)^) « 2{a^a) -|- I = 2n -|- 1. 
Therefore, the dashed curves in Fig. 4, representing the 
time evolution of the phonon number, also supply in¬ 
formation on the change of the oscillation amplitude of 
the resonator due to its coupling with the qubit. This 
phenomenon signifies the presence of a Majorana qubit. 
Certainly, for better performance of this hybrid system 
(e.g., with a higher fidelity), a higher quality factor of 
the resonator and a longer relaxation time of the qubit 
are preferred. 


VI. DISCUSSION 

Here we briefly compare our model to the one pro¬ 
posed by Kovalev et al. ,23 where a vibrating cantilever is 
utilized to rotate a Majorana qubit. The effective Hamil¬ 
tonian in their model [Eq. (7) in Ref. US] is in fact equiv¬ 
alent to the one in our manuscript [Eq. (5)]. This is 
understandable as both are in the framework of the spin- 
boson model. Note that for both cases there exists a 
static off-diagonal term [for our case, that is the 5 term 
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FIG. 4: (Color online) Time evolution of the occupations 
of the Majorana qnbit (solid curves) and phonons (dashed 
curves) which are in Rabi resonance. The qubit relaxation 
time is set as 50 ^s in (a) and 150 /rs in (b). The calcula¬ 
tions for both (a) and (b) are performed with two different 
resonator quality factors: Q — 10® and Q = 3 x 10®. 

in Eq. (5)] coupling the two levels of the qubit in the 
Hamiltonian. To neglect this term, in order to simplify 
the theoretical analysis, some conditions have to be satis¬ 
fied. Specifically, in Ref. 5, a certain equilibrium angle (0o 
there) of the vibrating cantilever has to be established. In 
our opinion, exactly solving this angle and then adjust¬ 
ing the experimental setup corresponding!}® are chal¬ 
lenging. However, in order to neglect the constant off- 
diagonal term in our case, the experimental setup must 
be mirror-symmetric about the middle point of the FM2 


gate, i.e., 1^2 = ^ 34 - Therefore, we think that our model 
is more easily accessible by experiments and hence more 
advantageous. 


VII. CONCLUSIONS 

In conclusion, we have proposed a hybrid system com¬ 
posed of a Majorana qubit and a mechanical resonator, 
implemented by a semiconductor nanowire in proxim¬ 
ity to an s-wave superconductor. In this proposal, three 
ferromagnetic gates are placed on top of and along the 
nanowire; the two outer gates are static and the inner 
one is free to oscillate harmonically as a mechanical res¬ 
onator. These ferromagnetic gates induce a local Zeeman 
splitting and give rise to four Majorana bound states, 
constituting a Majorana qubit in the nanowire. The 
dynamical hybridization of the Majorana bound states, 
arising from the motion of the oscillating gate, results in 
a coherent coupling between the Majorana qubit and the 
mechanical resonator. 

This hybrid system can be adjusted to be in resonance, 
e.g., with the assistance of a gate voltage on the ferro¬ 
magnetic gates, which controls the Rashba SOC locally in 
the nanowire. Our study reveals that under resonance, a 
strong coupling between the qubit and the resonator can 
be achieved. Consequently, quantum information can be 
effectively transferred from the Majorana qubit to the 
oscillator and then back to the qubit. This quantum in¬ 
formation transfer can manifest itself in modulating the 
motion of the oscillator, which may conversely signify the 
presence of the Majorana qubit. 
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